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We show that a hydrodynamical model with continuous particle emission instead of sudden freeze 
out can explain both the strange particle abundances and pion abundance from NA35 without extra 
assumption (e.g., sequential freeze out, modified equation of state, sudden plasma hadronization,...). 
In this scenario, the observation of a larger pion abundance is natural and does not imply a higher 
initial entropy and early plasma phase. 
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The main purpose of the ongoing and future heavy 
ion programs at the high energy laboratories CERN 
(Switzerland) and Brookhaven National Laboratory 
(U.S.A.) is to investigate the formation of hot dense mat- 
ter and the possible transition from hadronic matter to 
quark gluon plasma. Various possible signatures of the 
appearence of a quark gluon plasma (thereafter QGP) 
have been suggested: entropy increase (due to the release 
of new degrees of freedom, namely color), strangeness in- 
crease (due to enhanced strange quark production and 
faster equilibration), J/ip supression (due to color screen- 
ing), production of leptons and photons (emitted from a 
thermalized QGP and unaffected by strong interactions), 
etc. These signals have been studied extensively exper- 
imentally (see for exemple jjj). In this paper, we con- 
centrate on two signatures: strangeness increase which 
has been observed between p-p and nucleus- nucleus at a 
fixed energy and entropy production which is studied via 
pion measurements. 

A major problem to trace back any signature unam- 
biguously to a quark gluon phase is that it is still un- 
known which theoretical description describes best high 
energy nuclear collisions. On one extreme, one might use 
a microscopic model. Hadronic microscopic models fail to 
reproduce simultaneously strange and non-strange parti- 
cle data in nucleon-nucleon collisions and central nucleus- 
nucleus collisions at SPS energies (see ||,H ) ■ Partonic mi- 
croscopic models are expected to work at energies higher 
than SPS (however see Q and references therein). 

On the other extreme, one might use a thermal or hy- 
drodynamical model. In such models, it is assumed that a 
fireball (region filled with dense hadronic matter or QGP 
in local thermal and chemical equilibrium) is formed in 
a high energy heavy ion collision and evolves. Hydrody- 
namical models have been used successfully to describe 
various kinds of data at AGS and SPS. In particular they 
are able to account for strangeness data but in their sim- 
plest version, fail to predict big enough pion abundances. 

Since as already mentioned, both strangeness and pion 
or entropy productions are expected to be modified by 
the appearence of a quark gluon plasma, looking for a 
joint explaination (with or without plasma) of the rel- 
evant data is crucial. In the case of the thermal and 



hydrodynamical models mentioned above, various ways 
out of the pion problem have been proposed (see next 
section). In this paper, we study another possibility: use 
a more accurate emission mecanism, namely continuous 
emission, rather than standard freeze out, in a hydrody- 
namical description. This explaination has the advantage 
that no extra assumption is needed: once the initial con- 
ditions of the hydrodynamical expansion are fixed, both 
strangeness and pion yields come out with the right mag- 
nitude. 

Hydrodynamical or thermal description with 
(standard) freeze out emission - First we remind 
what is the status of the standard hydrodynamical or 
thermal description of relativistic nuclear collisions. In 
this kind of description, hadrons are kept in chemical or 
thermal equilibrium until some decoupling criterion has 
become satisfied. An exemple of freeze out criterion often 
used is that a certain temperature and baryonic potential 
have been reached. 

Since abundances are fixed by the chemical freeze out, 
the chemical freeze out parameters can be extracted by 
analyzing experimental particle abundances. This has 
been done by many groups ||. The models have some 
variations between them, but as a general rule, while they 
can reproduce strange particle abundances, they under- 
predict the pion abundance. This was first noted by |j] 
in a study of NA35 data and emphasized by J?],|| in an 
analysis of the WA85 strange particle ratios and EMU05 
specific net charge D q = (N + - N~)/(N + + N~) (with 
N + and N~, the positive and negative charge multiplic- 
ity respectively). 

Various possible mecanisms have been suggested so 
that a hadronic gas could yield both the correct strange 
particle ratios and pion multiplicity: sequential freeze 
out |9) or separate chemical and thermal freeze outs, 
hadronic equation of state with excluded volume cor- 
rections jn] |l2| , non-zero pion chemical potential , 
equilibrated plasma undergoing sudden hadronization 
and immediate decoupling Jq^|l3|]. 

All the physical points suggested to salvage the stan- 
dard hadronic gas model might need to be contemplated 
in a precise hydrodynamical description of relativistic nu- 
clear collisions, however it is somewhat suprising that one 
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has to go to such kind of details to reconcile the strange 
particle and pion data. Is it not possible to build a simple 
hydrodynamical model that yield the correct abundances 
without extra assumption? We discuss this question in 
the next section [ fL4| . 

Hydrodynamical description with continuous 
emission - In the standard hydrodynamical mod- 
els, one assumes that the freeze out occurs on a 
sharp three-dimensional surface (defined for example by 
T(x,y,z,t) = const). Before crossing it, particles have 
a hydrodynamical behavior, and after, they free-stream 
toward the detectors, keeping memory of the conditions 
(flow, temperature) of where and when they crossed the 
three dimensional surface. The Cooper-Frye formula fl5[ | 
gives the invariant momentum distribution in this case 



dN 



Ed 3 N/dp 3 



do-^fix.p). 



(1) 



da^ is a normal vector to the freeze out surface a and / 
the distribution function of the type of particles consid- 
ered. This is the formula implicitly used in all standard 
thermal and hydrodynamical model calculations of the 
previous section. 

The notion that particle emission does not necessarily 
occur on a three dimensional surface but may be contin- 
uous was incorporated in a hydrodynamical description 
in |[(|. In this model, the fluid is assumed to have two 
components, a free part plus an interacting part and its 
distribution function reads 



f{x,p) = ftree(x,p) + f in t{x,p). 



(2) 



/f rcc counts all the particles that last scattered earlier at 
some point and are at time a: in x. f- ln t describes all the 
particles that are still interacting (i.e., that will suffer 
collisions at time > x° and change momentum). The 
invariant momentum distribution is then 



Ed 3 N/dp 3 = d A xD^f tree (x,p)] 



(3) 



D /f ree (a: , p)\ is a covariant divergence in general co- 
ordinates and d A x is the invariant volume element. A 
priori formula ([|) is sensitive to the whole fluid history 
and not just to freeze out conditions as in formula (0). 

To compare particle abundances in the continuous 
emission and freeze out scenarios, we use a simplified 
framework to describe the fluid expansion, namely we 
suppose longitudinal expansion only and longitudinal 
boost invariance fTij] . This approximation allows to carry 
out some calculations analytically and turns the physics 
involved more transparent. It is implicit however that 
this description applies at best to the midrapidity re- 
gion. We will therefore concentrate on midrapidity abun- 
dances, precisely data from NA35. 

In this simplified framework, in the case of a fluid with 
freeze out at a constant temperature and chemical po- 
tential, the Cooper-Frye formul a (|lj ) can be rewritten ig- 
noring transverse expansion as |18| 
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(The plus sign corresponds to bosons and minus 
fermions.) It depends on the conditions at freeze out: 
and /if = [AbfoB + /j,s{ S, with B and S the baryon num- 
ber and strangeness of the hadron species considered, and 
MSfo^bfoi T{ ) obtained by imposing strangeness neutral- 
ity. So the experimental spectra of particles teach us in 
that case what the conditions were at freeze out. 

For continuous emission, we can approximate equation 
(0) as M 
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dyp±dp A 



y=0 ~ (2-) 2 

mi cosh rjTppdp 



d(f)dr] 



■ p± cos 4>pFTdr 



(5) 



/•p =0 .5 exp [(mj_ cosh 77 - ^i)/T] ± 1 

where V is the probabality to escape without collision cal- 
culated with a Glauber formula, tf (resp. pp) is solution 
ofP{T F ,p,(/),r);v ± ) = 0.5 (resp. V(t, p F , 0, -q; v±) = 0.5). 
In (||), various T and p = p^B + p$S appear (again p$ 
is obtained from strangeness neutrality), reflecting the 
whole fluid history, not just Tf D and pbfo- This history 
is known by solving the hydrodynamical equations of a 
hadronic gas with continuous emisison. It depends only 
on the initial conditions To and pbo- Therefore @ only 
depends on the initial conditions. 
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FIG. 1. Allowed region for the initial conditions deter- 
mined from the NA35 data. 

Once the spectra are known, they can be integrated 
to get abundances. Figure 1 shows the allowed re- 
gion of initial conditions that lead to the experimental 
NA35 midrapidity values @ A = 1.26 ± 0.22 Pi>0 . 5 GeV, 
A = 0.44 ± 0.16 Pi>0 . 5G ev, K° s = 1-30 ± 0.22 Pi>0 . 6 2GeV, 
h~ = 27±1 andp-p = 3.2±0.4. We do not use the K+ 
and K~ abundances because they were measured out- 
side the mid-rapidity region. For heavy particles, initial 
conditions dominate |2C| ] so when looking at their abun- 
dances in finite p± windows, transverse flow (assumed 
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zero initially) can be neglected as we did. For negatives, 
the experimental abundance is for all p± so tranverse 
flow does not affect their abundance either. The allowed 
window corresponds to T w 185 MeV, p,bo ~ 100 MeV 
, for an ideal gas equation of state and the strangeness 
saturation factor 7 S = 1.3 

Using a more sofisticated equation of state, the value 
of T might be decreased by some 10-15%, i.e., to 
155-165 MeV, compatible with (i.e., below) QCD lattice 
values for the phase transition temperature from QGP 
to hadronic matter. Our value of 7s is above 1 and this 
might looked surprising. However, its value is decreased 
by some 15% when looking at a more realistic equation of 
state. In addition, we have imposed strangeness neutral- 
ity, it is possible that this is a too strong constraint when 
analyzing data taken in a very restricted rapidity region 
(see plj where a similar problem was encountered). Using 
a larger value of 7 S , the size of the allowed window for 
initial conditions in figure 1 increases. There are other 
factors that influence the precise location and size of the 
window: values of the cross section needed to compute V 
(taken constant and equal for simplicity here), value of 
the cutoff in V = 0.5, etc. However the important point 
is that it is possible to find initial conditions of the hydro- 
dynamical expansion such that strange and non-strange 
particle abundances can be reproduced simultaneously 
without extra assumption. 

To illustrate why the continuous emission is able to 
reproduce both strangeness and pion data, in table 1, 
we compare results from the continuous emission sce- 
nario and the freeze out model. We took Tq = 185 MeV 
and /Zbo = 100 MeV for the continuous emission case, 
T fo = 185 MeV and fi bfo = 100 MeV for the freeze out 
case. We used 7 S = 1.3 for both p2fl . We expect roughly 
similar results for both models for heavy particles: in 
the continuous emission case, due to thermal suppres- 
sion, they are mostly emitted early [^pf , i.e., in similar 
conditions than in the freeze out model. Pions in the 
freeze out case are too few as discussed previously. Pions 
in the continuous emission case, on the other side, are 
emitted early and then on, so we expect to have more of 
them. This is precisely what we see in the table. 

The initial conditions that we discussed so far corre- 
spond to a hadron gas, starting its hydrodynamical evo- 
lution. However the present continuous emission scenario 
is in fact probably compatible with the possibility that a 
QGP was created before for the following reason. Con- 
tinuous emission is possible from a QGP but is inhibited 
by two factors: 1) only color singlet objects (not single 
quarks) can escape from it 2) the QGP core is surrounded 
by a dense hadronic region that the color singlets would 
have to cross without collisions to modify the previous re- 
sults. In this case, particle emission would occur mostly 
in the hadronic phase. Numerical estimates are however 
necessary to back up these qualitative arguments. 



t = 3fm 




T(MeV) 



FIG. 2. Massless pion gas: value of the pion density over 
entropy density as function of the initial temperature and po- 
sition for t = 3 fm in the continuous emission scenario (white 
surface). For comparison the value for the freeze out case is 
shown (gray surface). 

Relation between pion number and entropy - In 

a hydrodynamical model without shocks and dissipation, 
entropy is conserved. In the usual freeze out scenario, 
this is an important point because the initial entropy can 
be determined from the final multiplicity. To illustrate 
this connection, let us consider a massless pion gas. Sta- 
tistical mechanics yields the following relationship for the 
entropy density and the pion density s — 3.6/1^. In the 
longitudinal boost invariant model, the entropy conser- 
vation equation gives: st = const. Using dV — rdynR 2 
at y = 0, this can be written: dS{r) / (nR 2 dy) — const. 
So we can rewrite 



dS dS 
dy {T > ~ dy 



3.6- 



dy 



(6) 



Therefore, experimental knowledge of dN^/dy, permits 
to extract dS/dy, which are both independent of time. 

In the continuous emission case, s ^ 3.671^ because the 
distribution function needed to compute the pion density 
now depends on the escape probability V, embodied in 
/free (cf. eq. (2)). In fact, using similar methods than in 
|]l6| , one can show that 



3/3 -a 



4(1 + a) 



3.6/1^, 



(7) 



where a(t,p) = (4tt) 1 / d<j)d8 sinO [V/(l - V)] z=0 and 
f3(t,p) = (471-r 1 / d0d9 sin 9 cos 2 9 [V/(1-V)] Z=Q . This 
relationship is plotted in figure 2. For all radii p, time r 
and inicial condition To, s(r, p) < 3.6n^(T, p). One then 
gets in the longitudinal boost invariant model: 



dS <Z6 dN„(r) 
dy dy 



(8) 



in other words, for a fixed entropy, there are more pions 
in the continuous emission case than in the freeze out case 
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(for all times t). This is in fact expected (even using a 



mor p realistic equation of state and including pions from 
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decays): since pions are emitted continuously, there are 
more copious than in the the usual freeze out case for a 
given dS/dy (itself fixed by the initial conditions). 

As a consequence, a large experimental value of 
dN^/dy should not necessarily be associated to a large 
dS/dy and be considered a hint of QGP formation as 
usually done in freeze out models M. A larger dN^/dy 
and not large dS/dy is a natural consequence of contin- 
uous emission, compared to freeze out. This result is in 
contrast for exemple with JtLH. 

Conclusion - In this paper, we discussed data on 
strange and non-strange particles by NA35, from an hy- 
drodynamical point of view. The standard model with 
sudden freeze out can reproduce the strange particle data 
but underpredicts the pion abundance, if no extra as- 
sumption is made. We showed that a hydrodynamical 
model with a more precise emission process, continu- 
ous emission, can reproduce both the strange and non- 
strange particle data without extra assumption. 

This indicates the necessity of doing a more accurate 
description of particle emission in hydrodynamics, else 
some problems might artificially appear, such as a too 
low predicted pion abundance as discussed here or a too 
high freeze out density 0] . 

This point is reinforced by the fact that a large pion 
number is usually associated with a large entropy and 
QGP formation. Here we showed that a large pion num- 
ber can be generated by continuous emission without 
modifying the entropy. (The larger pion emission will 
cause a faster cooling and shorter fluid lifetime.) In other 
words, a better understanding of particle emission in the 
hydrodynamical regime is also necessary to assess the 
possibility of QGP formation in relativistic heavy ion col- 
lisions. 

There is a growing tendency to use hydrodynamical 
models to describe relativistic nuclear collisions, there 
is also a growing concern to modelize freeze out better 
p3| , ^4| . However it seems that the very notion of particle 
emission during the hydrodynamical expansion needs to 
be put under more scrutiny. 
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TABLE I. Comparison of experimental particle abun- 
dances with continuous emission and freeze out predictions. 

Experimental value Continuous emission Freeze out 
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